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Abstract 



Recently the scaling function of the dilute non-contractible self-avoiding 2D polymer 
loop on a cylinder was related to the Painleve III transcendent. Using the perturbation 
J> ■ theory, the thermidynamic Bethe ansatz and numerical calculations we argue a similar 
relation for the contractible self-avoiding loop. 



1. Introduction 



Among the recent studies in the 2D integrable relativistic field theory (RFT) one of the 
most impressive results was achieved by P.Fendley and H.Saleur in ref. [1] (see also related 
^ . ref. [2]). The authors have succeeded to link the universal scaling funcition of a single 2D 
self-avoiding loop winding once around a cylinder to a particular Painleve III transcendent. 
I After almost 20 years of the famous two-spin Ising scaling correlation function [3] standing 
^ ' alone, a non-linear differential equation appears again in 2D RFT to govern the universal 
scaling behavior. Being related to an apparently interacting field theory the result of 
[1] seems extremely suggestive. Moreover, it solves exactly an important problem in 2D 
polymer statistics and is directly comparable with series expansions, simulations etc. It 
should be mentioned that the discovery was preceded by a series of fascinating studies in 
the topological and N = 2 supersymmetric 2D RFT's [4-6]. Although it is not obvious 
that it is the N = 2 SUSY that plays the most important role in the phenomenon of [1] , the 
physical significance of ref. [4] worth to be better understood (see ref. [2] in this connection) . 

The problem of 2D self-avoiding polymers is formulated basically as follows. Consider 
some 2D lattice (say the honeycomb one to get rid of any problems with self-avoiding) and 
count the continuous self-avoiding paths (closed or open) through the links of the lattice. In 
general there are several disconnected components, each of them being called the polymer 
(respectively closed and open). Because of the self-avoiding there are many topologically 
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different configurations and tlierefore an innumerable amount of various counting problems. 

As a first step one tries to count the closed polymer configurations taking into account 
the number of separate polymers. The configurations can be generated e.g. by the K- 
expansion of the 0(n)-model partition function 

Z{K,n) = ir J] + (1.1) 

<i,3> 

where the product is over the nearest neighbour pairs of sites < i,j >. The n-component 
spins s"; a = 1, 2, . . . , n sit at the sites and are chosen so that trs'* = 0, trs^s** = 5°"'^ and 
trs^s^'s'^ = 0. For (1.1) we have 

Z{K,n)= K*°^ ^'''^^n*°^ (1.2) 

loop 
COnfiguraitons 

As usual for a homogeneous lattice of large volume V the thermodynamic limit is supposed 
to exist 

J^{K,n) = -^ lim log Z{K,n) (1.3) 

defining the specific free energy J-'{K, n). 

What is interesting for the field theory is the statistics of large loop configurations 
(i.e. with the number of lattice links in a polymer tending to infinity) which is believed 
to be universal, i.e. independent on the lattice realization. The loops blow up near some 
critical value Kc where the statistical quantities (like the free energy) develop singularities 
characteristic for the second order phase transition. The phase at K < is called the 
dilute polymer one. While K < the correlation length (i.e. the typical spatial extension 
of the polymers), which we denote here, remains finite but diverges as 

~ {Kc - (1.4) 

The standard parameter p which appears in the exponent of eq.(1.4) is related to n as 
follows 

TT 

n — 2 cos — (1-5) 
P 

Near the critical point the observables (correlation functions etc.) apportion some singular 
in Kc — K parts, which depend only on the distances and extents scaled in the units of 
and bear the rotational (Lorentz) symmetry characteristic for RFT. Being universal 
these scaling functions are objects of the field theory. E.g. the free energy (1.3) itself 
behaves as 

T{K,n) = T{K,n) + M'^f{n) (1.6) 

where is some non-universal background regular m. K — Kc which is typically of no 
interest for the field theory. Contrary, the second piece in (1.6) is contributed by the 
asymptotically large polymer loops and describes their universal statistics independent on 
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how the polymers are arranged on the lattice level. For the 2D polymer problem f{n) is 
known exactly 

/(n) = -ltan^ (1.7) 

This result comes both from the relation to the sin-Gordon model (see below) where the 
exact vacuum energy was found by different methods [7-9] and from a formal analytic 
continuation of the thermodynamic Bethe ansatz (TBA) equations at integer p [10]. 

From eq.(1.2) it is clear that the n-expansion 

/(„) = ^ = -ln-ln^-^(l+'jf)n^ + ... (1.8) 

k 

is in fact the expansion in the number of separate polymer loops. E.g. the first number 
/i = —1/4 corrsponds to the (scaling) internal statistical weight (or activity) of a single 
isolated loop. In other words /i is the specific free energy of a dilute gas of non-interacting 
contractible polymer loops. The second two-polymer number /2 = — l/47r receives in fact 
contributions from two topologically different polymer configurations drawn in fig.l. The 
first one (fig. la) gives the first virial coefficient for the gas of contractible loops which takes 
into account their interaction (due to the self-avoiding) at the two-component cluster level. 
The second contribution (fig. lb) is the activity of the structure with two nested loops. 
Further numbers in (1.8) include more and more distinct topological contributions. It is 
an interesting problem (unsolved to my knowledge) to separate them. 

However, in the present note we deal with another settlement of the polymer problem. 
In the next section the counting problem of closed polymers on an infinite cylinder is 
considered and the corresponding scaling functions are defined. Then the remarkable result 
of ref.[l] is quoted which relates explicitly the one non-contractible polymer scaling function 
to a particular Painleve III transcendent. We present also a similar relation for the scaling 
function of a single contractible polymer. Sect. 3 contains few details and propositions about 
the Painleve III transcendents including their representations as Fredholm determinants 
as well as the TBA- like representations [1,2]. In sect. 4 we discuss the useful relation 
between the cylinder polymer scaling functions and the finite-temperature free energy of 
the sin-Gordon model (SG) in the repulsive regime > Att. In particular this includes 
the restricted sin-Gordon models (RSG), i.e., the $13 perturbed minimal conformal field 
theories (CFT's). Few first terms in the perturbative expansions of the SG and RSG 
finite-temperature free energy are developped in sect. 5. They are to be compared with 
the corresponding expansion of the Painleve III. Sect. 6 contains the TBA considerations 
in the sin-Gordon model near its = 2 supersymmetric point = 167r/3. These permit 
us to derive the TBA-like formulas quoted in sect. 3. A bulk of questions and hints arises 
in connection with the things considered. Some of them are listed in sect. 7. 



2. Polymers on a cylinder. The Painleve III 

Let us start again from the lattice level. Imagine a long cylinder made of say again 
the honeycomb lattice (fig. 2) and denote L the length of the cylinder and R its circum- 
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ference. When counting the closed polymer configurations on this cylinder lattice one can 
distinguish between two kinds of polymers, the "winding loops" which wind once around 
the cylinder (fig. 3a) (it is plain that the self-avoiding allows a closed polymer to do it only 
once) and the "non- winding loops" which don't (fig. 3b). We shall count these two kinds 
separately defining the partition function as 

Z{m, n\K,R,L)= ^ K*°^ lmks^#of winding loops^#of non-winding loops (2.I) 

cylinder polymer 
configurations 

where m is the weight of a winding polymer and n is that of a non-winding one. At L — > oo 
as usual we introduce the free energy per unit length 

£{m,n\K,R) = lim log Z{m,n\K,R,L) (2.2) 

Of course, at R fixed (in the lattice units) £{m,n\K, R) shows no criticality near Kc- 
It does however if R simultaneously goes to infinity as fast as the correlation length 
does. Defining the scaling circumference 

t = MR (2.3) 

we expect that at t fixed and K ^ Kc — 

S{m, n\K, R) = S{m, n\K, R) + MF{m, n\t) (2.4) 

Here £{m, n\K, R) is non-singular at Kc and in fact for R large enough 

£{m, n\K, R) = RT{K, n) (2.5) 

independently on m (with the same T{K,n) as in eq.(1.6)). What is interesting for us 
is the universal scaling function F{m,n\t) which describes the statistics of large (dilute) 
polymers on the cylinder. Obviously at t — > oo (i? » M~^) the influence of the cylinder 
geometry disappears and 

F{m,n\t) ^ tf{n) , t ^ oo (2.6) 

with /(n) defined in the previous section. For what follows it is convenient to introduce 
also the special notations 

F+{n\t) = F{n,n\t) (2.7) 
for the scaling function of equally weighted winding and non-winding loops and also 

F_{n\t) ^ F{-n,n\t) (2.8) 

for that weighted with the opposite signs. 
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Again the double expansion in m and n 

F{m,n\t) ^J2^ii,i2iiW'n^' 

h,h (2.9) 

= mFi,o(t) + nFo,i(t) + m^F2,Q{t) + mnFi,i(t) + n^Fo,2(t) + • • • 

isolates the configurations with fixed numbers li of winding loops and I2 of non-winding 
ones. Each of the linear terms i^i,o(^) and -Fo,i(0 corresponds respectively to the scaling 
statistical weight (per unit length of the cylinder) of a single isolated non-contractible 
(fig.Sa) and contractible (fig.Sb) polymer loop. The next quadratic terms in expansion 
(2.9) contain different two-loop scaling functions. Namely, -F2,o(^) and i^i,i(t) are related 
to the two-polymer clusters drawn in figs. 4a and 4b while i*b,2(^) is contributed by two 
topologically different configurations of figs. 4c and 4d. 

The one-loop cylinder scaling functions are of primary interest below so that we in- 
troduce special notaitons for them 

Fjt)=Fio{t) 

In connection with the functions (2.7) and (2.8) we shall also consider the combinations 

F±{t) = F,{t) ± Fn{t) (2.11) 
The remarkable result of ref.[l] is that 

F„W = ^ (2.12) 

where U {t) is a particular solution (regular at t > 0) to the Painleve III equation 

1 d d 1 

-—t—U=-smh2U (2.13) 
tdt dt 2 ^ ^ 

This equation admits a one parameter family of regular at t > solutions (see e.g. [11]) 
called the Painleve III transcendents. The special solution in eq.(2.12) is fixed by the 
following boundary conditions at t — > 

C/(t) = -^logt-ilog^ + 0(t^/3) (2.14) 

where we have denoted 

«=P2|^ = 3.91392514... (2.15) 

Prom eqs.(2.13) and (2.14) it is quite easy to see that U{t) is essentially a regular series 

1 1 °° 

U{t) = -4 loge - 2 log ^ + J] UnC (2.16) 

n=l 
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in the variable 

e = t^/^ (2.17) 
The coefficients t/„ can be recursively generated by the equation (2.13) 
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f^nV.r = ^e^Sl.^"^" - ^e-^E:..^^r ^2.18) 



9 ^ K 16 

n=l 

This way one finds 

1 4 °° 
tF„(t) = -- + -^nt/nr 

n=l 

_ 1 3 3(^«3-i8) 2 27(/t3 + 18) 3 27(av6 + 30k3-648) 4 
~~3^4^^ 128^2 ^ 2048^«3 131072^^ ^ 

(2.19) 

On the other hand the Painleve III transcendent U {t) admits the following infinite series 
representation [11,2] 



oo „2fc+l -t cosh Si JQ 



where it is implied that 92k+2 = ^i- The series is convergent for all t > while the first 
terms control the t ^ oo asymptotics. In particular the scaling function F„(t) falls off 
exponentially at t — > oo 

Fr,{t) ^ ~K^{t) + O {e-'') (2.21) 

where 

K^(t) ^^j e'^^-'^^'^'^de (2.22) 

is the modified Bessel function. The non-contractible one- loop function F„(t) is plotted in 
fig.5. 

Below we argue that the contractible one-loop scaling function Fc{t) satisfies the 
following equation 

-^tFc{t) = -- cosh2C/ (2.23) 
where U{t) is the same Painleve III transcendent (2.20) as for the non-contractible function 

Fn{t). 

Since U{t) ^ at t — > oo it follows from eq.(2.23) that Fc{t) has a linear leading 
asymptotic in this limit 

Fc{t) ^ at t ^ oo (2.24) 
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in agreement with (2.6) and the expansion (1.8) of the bulk free energy. Eq.(2.23) defines 

Fc{t) up to an item of C/t with an arbitrary constant C. It seems natural to choose C so 
that the corrections to (2.24) at t — > oo would be exponentially small. Then we have 

Fc(t) = -^ + ^y^ K'o{r)rdr + O {e-^') (2.25) 

In the next section it is argued that in fact Fc{t) develops a series representation quite 
similar to eqs.(2.12) and (2.20). Define V{t) as 

2 

where again ^2fc+i = Oi in each term in the r.h.s. (note that now we sum up the "even 
loops" instead of the "odd" ones in eq.(2.20)). Then 

F.(t) = -"^ (2.27) 

At small t the contractible loop scaling function is again a regular series in the variable 
(2.17) 

4 °° 

tF,{t) = c--Y,nVnC (2.28) 

n=l 

The coefficients Vn are readily restored from eq.(2.23) 

9 ^ ^ K 16 ^ ' 

n=l 

As it was mentioned before, the constant term c in (2.28) in not fixed by eq.(2.23). We shall 
see shortly that the above choice of the integration constant C in eq.(2.25) corresponds to 
c = 1/18. Expansion (2.28) reads 

1 3 3(^3 ^ 18) , 9(5^3 _ 'XJi^f + Q^ + m^ 
^^'^^^^ - 18 ~ 4^^^ ~ 128/.2 ^ + 2048/.3 ^ 131072;^^ ^ ^^"^^^ 

The contractible scaling function Fc{t) is presented in fig.6 (without the "bulk energy" 
term -t/A). 

For the later use we quote also the short-distance expansions of the combinations 



and 



^ ^^ 7 3, 27 ,o 9(^3- 54) ,o 81(^3 -54),. 
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These two functions are plotted in fig. 7. They satisfy the equations 

tdt ' 2 



(2.33) 



3. Painleve III vs. Fredholm theory 

We begin with the second kind Fredholm equation 

/oo 
K{e,e')f{e')de' = g{e) (3.1) 
-00 

where K{9,9') — K{9',9) is a symmetric kernel (bounded in L2(— 00,00)) and A is the 
spectral parameter. Below we take the kernel of the following special form 

-uie)-u{e') 

mo') = - — jzrr (3-2) 

cosh 

2 

with some suitably chosen "potential" u{9). For the Painleve III theory 2u{9) — tcosh.9 
will be relevant. Note, that in terms of x = and substituting f{9) — exp{9/2)f{x), 
g{9) = ex.p{9/2)g{x) eq.(3.1) reads 

/■oo —u(x) — u(x') 

fix) - 2A / ^^f(x')dx' = ~g{x) (3.3) 

Jo x + x' 

The spectral information about (3.1) is contained in the Fredholm determinant -0(^1 A), 
which is entire function of A with zeroes located at the eigenvalues A^, a = 1, 2, ... of (3.1). 
It can be calculated as a convergent series in A 

= E M^^n(«) (3.4) 



n=0 



with 



DJu) = / det f ^ f] 2e-^''^''^Ux 

Jo \Xi + XjJf-J: 



I n 



i=l 

2 



(Xi Xj) YJ ^-2u(xi)^^i 
( 



(3.5) 



TT^-2«(a;i)^ 
{Xi + X^Yl}^ 



e 

Xi 



Denote 



-W{u\\) = log D{u\X) (3.6) 
For M^(u|A) we have the following A-series 

W{u\X)^Y.—^r.{u) (3.7) 

n=l 



convergent while |A -"^l > sup^ \^a'^\- Here 

-2u{ei) POO n r)-2u{xi) 

g-^de. - / n 



Wn{u) = / n 3 n = / n -— dXi (3.8) 



where we imply that 9n+i = Oi or a^^+i = ^^i- 

A link to integrable non-linear equations appears if we take 



m = le'+le-'+m (3.9) 



with some r{9) and denote 



U{u\X) = W{u\\) - W{u\ - A) = log (3.10) 

Then it happens that (we denote d — d/dt and d = d/dt) 

4.ddW{u\X) = ^ (e^u(^\^^ - l) (3.11) 

or in terms of the determinants D[u\X) 

D{u\X)ddD{u\X) - dD{u\X)dD{u\X) = ^ {D^{u\X) - D^{u\ - A)) (3.12) 

In particular U = U{u\X) = —U{u\ — A) solves the 2D sinh-Gordon equation 

4aaC/= ^sinh2t/ (3.13) 

Eq.(3.11) can be verified comparing order by order the A-expansions in the right and left 
hand sides. This would require the following relations between the integrals (3.8) 

15agw„w= E n^^T'^^T^' (3.14) 

f Pl,P3.P5.->0 \k=0^ '^^ ^ ^ 
Ipi+3P3+5P5 + ...=«/ 

E.g., 

^^^Wx = Wi 

AddW2 = 

9 



etc. The first three relations in (3.15) are verified directly by comparing the corresponding 
integrands. Starting from n = 4 one has to perform a symmetrisation in the rational in 
Xi, i = 1,2, ... ,n part of the integrands (this is allowed due to the permutation symmetry 
of the rest part J|"^]^ exp(— 2u(xi))(ia;i, common for the both sides of eq.(3.14)). I have 
checked (3.14) up to n = 6. The complete proof perhaps can be found following the lines of 
ref.[ll]. It should be stressed that the validity of (3.14) is of purely combinatorial nature, 
being based on certain algebraic identities between the symmetrized rational expressions 
in Xi. Therefore we expect eqs.(3.11-14) to hold independently on A and even on the choice 
of the "residual" potential r{9) in eq.(3.9). Moreover, they remain valid if in eq.(3.1) (and 
therefore also in eq.(3.8)) we replace the whole real axis by any other sensible integration 
contour in ^. In this note we shall not develop further these lines of generalizations. 

Before turning at the rotationally symmetric (isotropic) version of eq.(3.11) it is worth 
to mention the following interesting observation of refs.[l,2]. Consider a small variation 
6u{9) of the potential in (3.2). Then 

SWiu\X) = -2A J R{e\X)Su{9)d9 (3.16) 

where R{9\X) = 71(9, 9\X) and n{9,9'\X) = n{9',9\X) is the resolvent kernel of eq.(3.1), 
i.e., the (unique at A 7^ Ao) solution to 

n{9, 9'\X) -X j K{9, 9")n{9", 9'\X)d9" = K{9, 9') (3.17) 



As a series in A it reads 



R{9\X) = ^X''Rn{9) (3.18) 



n=0 



where 



Rn{0) = e-^<'^ / 3 2 3 ...d9r^ (3.19) 



-2u{0i)-...-2u{0r,) 

cosh - — r— ^ cosh , . . . cosh 

We define also 



2i?+(^|A) = R{9\X) + R{9\ -X) = 2J2 >^^^R2k{0) 

k=0 

oo 

2R-{9\X) = R{9\X) - R{9\ - X) ^ 2Y,X^^^^R2k+i{9) 



(3.20) 



k=0 



In refs.[l,2] on the basis of the field theory considerations quantity R^{9\X) was related to 
the following "TBA-like" system of non-linear integral equations 



10 



with the same potential u{9) and the same A. It was argued that there is a solution e{9), 
r]{e) to (3.21) such that 

R+{e\X) = e-^W (3.22) 



In sect. 6 of this note we extend slightly the TBA considerations of [1,2] to argue a similar 
formula for R_{9\X) 



Strictly speaking, the arguments of [1,2] (and also of sect. 6 below) are relevant only 
for the isotropic case 2u{6) = tcosh^. However, if we start to solve the system (3.21) as 
a power series in A we would encounter for the corresponding terms in (3.22) and (3.23) 
expressions like 



Here the "potential weighted" integrations over 9i,...,9n arise from the iterations of 
eq. (3.21b) while /(6'i, . . . , 6^) results from the "intermediate" integrations in the first equa- 
tion (3.21a) or in eq.(3.23). They turn to be rational functions of Xi = exp Of, i — 1,2, ... ,n 
which, after symmetrization, are expected to reproduce the (symmetrized) rational parts 
in the integrands of (3.19). In first few orders this can be verified explicitly. If (3.21-23) 
is true in the isotropic case (as it was argued in [1,2] and sect. 6) this phenomenon must 
persist at all orders, since it is quite hard to imagine that the transcendental integrations 
over 01,..., On may play any role in relating (3.24) to (3.19). One concludes that the 
TBA representation (3.21-23) is of more general nature and must hold for any reasonable 
potential u{9). Also, corresponding to different possible choices of the integration contour 
in (3.1) we have to change the contour in (3.21b). Since the ^-integrations in eqs. (3.21a) 
and (3.23) are important in the building of I{9i, . . . , 9^), the contours there must remain 
unchanged. Of course, the consideration above is not a proof, which is lacking up to now. I 
have verified the TBA representation numerically for several (presumably random) samples 



The TBA representation (although unproven) turns extremely useful in the numerical 
calculations. E.g., for the isotropic case (3.25) with |A| < 27r (and not too close to 2tv) the 
iterative solution of the system (3.21-23) converges fast and works much better then both 
the Painleve equation and the direct calculation of the Fredholm determinant. This typi- 
cally persists for other samples of u{0) with |A~-^| well below the sup^ \^a'^\- In particular, 
the curves of figs. 5-7 are computed in very this way. 

From now on we switch to the rotationally symmetric (i.e., invariant under t — > At, 
t A~^t) solutions to (3.11-13). This corresponds to taking in (3.9) r{9) = const (which 
can be absorbed by the spectral parameter) and keeping (— oo, oo) as the integration region 
in (3.1). Thus we restrict to t = t real and 




(3.23) 




(3.24) 



ofu{9). 



2u{9) = t cosh 9 



(3.25) 
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In eqs.(3.11-13) we can substitute 2d=2d = d/dt, Add = d^ /dt^ +t-'^d/dt so that (3.13) 
is reduced to the Painleve III equation (2.13). The series (3.7) turns also the large t 
expansion and we have 



W{t\X) = 2XKo{t)+8X^ 



I'OO I'OO 

\ogt I KQ{r)rdr + / KQ(r)rlogrdr 
Jt Jt 



+0{X^e-^^) (3.26) 



Considering real solutions we have to take A real. It is possible to show that at real t > 
sup„A-^ < 27r. Therefore W{t\X) is regular at t > for A < (27r)-^ At 27rA > 1 the 
Fredholm determinant D{t\X) has infinite number of zeroes on the positive t axis, t = 
being their condensation point [11]. We shall consider the case 27r|A| < 1 where all the 
functions appearing in eqs. (3.6-13) are regular at t > 0. Denote 



^ . . TTCr 

27rA = sm — — 
2 



(3.27) 



with real — 1 < cr < 1. It is relatively easy to analyse the t — > behavior of the determinant 
(3.4) with the potential (3.25). One finds 



W{t\X) = log ^ + B{a) + O {t'^'-) 



(3.28) 



where {^^{x) — r'{x)/r{x)) 



(l+x) 



1 + x 



-1-x 



- 2 



dx 



(3.29) 



In particular, for the combination (3.10) we have 



1 



U{t\X) = a log - - - \ogK{a) + O {t^^^'^) 



(3.30) 



with 



K{a) = exp(2S(-(7) - 2B{a)) = 



1 a 
2+2 



(3.31) 



Using (3.30) and eqs. (3. 11-13) one can systematically recover the further t — > corrections 
to (3.28). They appear as 



W(t|A)= <"+f -' logg-g+ f; B„..„( 

m,n=0 

where am,n — {m + n) + a{m — n) . In this double series 

Bo,o{(t) = B{a) 
12 



2ar 



(3.32) 



(3.33) 



and the next coefficients 5^ „(cr) are found order by order from the following relation 



In 



,2ao,i 



m,n=0 



k((7) 



exp 



2 ^ [-Bm.nlo") - -Bn,m(-0')] 



X 



2a, 



v 



m,n=0 
(m,ri)7^{0,0) 



Obviously Bm,oicr) = for m > 0. Few first Bm,nS are 

5o,i(^) = 



So,2(fT) = 



/.(cT)(l-a)2 

8 

8 

«2((7)(l-(7)4 



(3.34) 



(3.35) 



1/3. 



The polymer one-loop scaling functions of sect. 2 correspond to the special case a = 



4. Relation to sin-Gordon 

The sin-Gordon model is a theory of 2D scalar field (fi{x), = (x^,x^), with the 
action 

AsG = ^ J{d^(pf(fx -fx J : cos Pep: (fx (4.1) 

where /3 is a real dimensionless parameter (/3^ < 87r). Here : . . . : denotes the normal 
ordering with respect to massless free fields. To be precise we choose such an infrared 
cutoff that at = the two-point correlation function < ip{x)Lp{0) > is 



(V?(a;)(/p(0))^=o = -^log|a;| (4.2) 



In this scheme the (real) coupling constant /i is dimensional /i ~ [mass] with 
Sometimes it is convenient to use another (positive) parameter 
instead of 
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The action (4.1) is obviously symmetric under the global shifts of the field (p 



if{x) ip{x) + 



27rn 



(4.5) 



where n is an a:-independent integer. Therefore the states > of the model can be 
classified in their behavior under (4.5). The irreducible states |\E'q, > are characterized by 
the index (quasimomentum) —13/2 < a < 13/2 such that 



2ti 



(4.6) 



For each a there is an infinite dimensional space of states Aa-, the complete space being 



A= 

-l3/2<a<(3/2 



(4.7) 



In the infinite volume system the transitions (4.5) are suppressed to zero so that all the 
spaces Aa are completely degenerate in a. However in a finite geometry (we imply that 
the finite geometry settlement respects the symmetry (4.5)) these transitions are allowed 
and a becomes an important parameter of the theory. Let us call Aa the a-sector and the 
corresponding ground state \^a^ > the a-ground state. The a-ground state energy Ea is 
a periodic function of the quasimomentum a ^ a + For symmetry reasons we expect 
Eoi to have a minimum at o; = (mod /3Z), i.e. I^q'^ > is the true ground state. 

Before considering the finite geometry effects remind briefiy the well known struc- 
ture of the infinite volume space A^ which is essentially independent on a. Aa contains 
the infinite volume ground state (the ct-vacuum) together with the excitations which are 
(massive at jJL ^ in (4.1)) relativistic particles subject a to factorized scattering. The 
spectrum of particles always contains the soliton-antisoliton doublet (s,s). Its mass M is 
related to the coupling ji in eq.(4.1) as follows [8] 



P 



p + l 



p+ 1 



VttMF 



p+l 



-, 2/(P+l) 



2r 



(4.8) 



Quasiclassically one can think about s and s as of the kink configurations of field ip{x) 
such that (currently we imply to be the spatial coordinate) 



(p{x^ oo) — (f{x^ — > — oo) = ±27r//? 



(4.9) 



for the soliton and antisoliton respectively. At p < 1 the spectrum includes also the s 
bound states n = 1, 2, . . . < 1/p of masses 



rrin = 2M sin 



TTpn 



(4.10) 
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The factorized scattering amplitudes of all these particles can be found e.g. in [12]. Note 
that all the observables are independent on a. The (specific) ground state energy [7,8] 



£: = -^tan^ (4.11) 
4 2 ^ ^ 

is also a independent. 

Let us now put the sin- Gordon model on a finite space circle of circumference R and 
impose the periodic boundary conditions, i.e., (^{x*^, + R) — ip{x'^,x^). In the euclidean 
version of (4.1) this corresponds to the geometry of infinite (or very long L — > oo in the 
"time" direction x^) flat cylinder. We denote Ect{R) the corresponding ct-ground state 
energy. In the perturbed conformal field theory (CFT) picture this observable would 
correspond to the specific (per unit length of the cylinder) free energy with the scalar 
operator exp(iQ;(/?) "flowing along the cylinder". The UV conformal dimension of this 
operator is Aq, = ck^/Stt. Therefore we expect at — > [13] 

In general 

E^{R) = MFsG{a,p\MR) (4.13) 

where the soliton mass M is used as an overall scale to isolate the scale independent function 
FsG and the dependence on the sin-Gordon parameter (4.4) is explicitly indicated. From 
(4.12) we have at t ^ 

FsG{a,p\t) ^-^Jl-^], t^O (4.14) 



6t \ 7T ^ 

independently on p. In the opposite limit t — > oo we expect from (4.11) 

t TVp 

FsG{c(,p\t) — — - tan — -|- corrections, t — >• oo (4-15) 

independently on a. The leading correction in (4.15) comes (at p > 1/3) from the virtual 
soliton or antisoliton trajectory winding once around the cylinder. In view of (4.9) it is 
plain that in the ct-sector these trajectories are weighted by the factors of exp{±2i7ra/ P) 
respectively. Therefore, summing up the s and s contributions 

t np 2 27ra , , , 
FsG{ct:,p\t) + - tan — = cos — -— ivi(t) -|- (multiparticle or bound state contributions) 

4 2 TT p 

(4.16) 

Considering at the microscopic level one can argue [14] that the described above 
cylinder sin-Gordon settlement is in the same universality class as the cylinder polymer 
counting problem of sect. 2 provided the non- winding polymer weight n is related to the 
sin-Gordon parameter (4.4) by eq.(1.5). Also, comparing (4.16) with what one expects at 
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large t for the polymer counting scaling function F{m, n\t) we can relate the sin-Gordon 
quasimomentum 2TTa/ (5 to the winding polymer weight m 



27rQ; / , , „x 

m = 2cos— - (4.17) 

In addition it is clear that the soliton mass M has to be identified with the inverse cor- 
relation length of the polymer problem. Altogether it seems not misleading to think of 
the polymer loops as of the virtual soliton trajectories (summed over the "orientations" 
s and s). Note, that the sin-Gordon-polymer relation (1.5) implies that p > 1 where the 
sin-Gordon spectrum contains no extra particles but s and s. It looks quite natural that 
the self-avoiding polymers form no bound states. 

We conclude that if m and n are related to p and a as in eqs.(1.5) and (4.17) 

F(m,n|t) = F,G(a,p|t) (4.18) 

This correspondence turns rather useful. It permits often to reinterpret the sin-Gordon 
results in the polymer terms. E.g., eq.(1.7) is simply read off from the ground state energy 
(4.11). Vise versa (4.18) implies that at m = n = 0, i.e. p = 2 and 2a/ (3 = 1/2 the 
sin-Gordon a-ground state energy FsGici,p\t) vanishes for all t. This is indeed the case 
due to the N = 2 supersymmetry of the sin-Gordon model at this point [15]. One then can 
expand FaQ{a,p\t) around this point to obtain the scaling functions Fi^^i^{t) in eq.(2.9). 

The standard tools of studying the sin-Gordon model become applicable to the poly- 
mer problem of sect. 2. In particular, in the next section we use the perturbation theory of 
(4.1) in n to evaluate the small t behavior of F(m, n\t). The perturbative series for Ea{R) 
is constructed in the usual way. In fact it is the expansion 

oo 

^a(i?) = -^E^2n(a)r7^^ (4.19) 



6R 

n=0 



in the dimensionless variable 



r] = -27r// ( — ) = -27rrp ( — ) (4.20) 



^ \ 2/(p+l) / ^ \ 2/(p+l) 

where 

rp = //M-2/(f+i) (4.21) 
is the "//-M ratio" defined by eq.(4.8). The first coefficient in (4.19) 

co(a) = l (4.22) 

TT 

is the efi'ective central charge corresponding to operator exp{ia(fi) while the subsequent 
ones are given by perturbative integrals. Introducing the complex coordinate z = + ix^ 
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on the cylinder 2; ~ 2; + i?2 and substituting u = exp{27riz / R) we have explicitly 

12 

C2n(Q;) 



(«i = l) 



(27r)2"-i(n!)2 

n 



X 



\Ui — Un 



Vi - V, 



|4A 



n 



|4A 



n 

i=l 



. |2A(l-2a//3) I |2A(l+2a//3) 



subtrac- 
tions 



(4.23) 

where the disconnected parts are subtracted to build the connected 2n-point cylinder 
correlation function. The integrand is invariant under the homogeneous shifts along the 
cylinder Ui Aui, Vi Avi so that we can fix one of the points {ui in eq.(4.23)). Note 
that the perturbative coefficients C2n(ct) are not periodic in the quasimomentum a. In fact 
expansion (4.19) determines the a-ground state energy only if \a\ < /3/2. For \a\ > P/2 
and a = ao (mod /3), |q;o| < /?/2 it corresponds instead to an excited state in the ckq 
sector. The discontinuities occur at a = nP/2 with n integer, where the unperturbed 
space of states contains two (connected perturbatively) degenerate states exp(±m/3(/?/2) 
and expansion (4.19) has to be modified. 

The series (4.19) admits another interpretation if 2a/ (3 — ±(l/p — k) with k a non- 
negative integer [16] (note that in this series only k = and 1 would correspond to a-ground 
states). Take e.g., 2a/ (3 = 1/p — k and consider first the integration over Vj, z = 1, 2, . . . , n 
in eq.(4.23). Following [17] we have 



1 



n 



Ui - u. 



|4p/(p+l) 



Wi - V, 



|4p/(p+l) 



.|4p/(p+l) 



J] |„.|2(fcp-2)/(p+l)|^.|-2/cp/(p+l)^2^. ^ 



i=l 



(4.24) 



2n 

Qp 



($l,fe+l(0)$l,3(«l)^l,3(«2) . ..^l,3M^hk+l{oc))M, n l^il"^^^''^' 



where $ifc,/c = l,2,... denote the (thermal series) primary fields in the minimal CFT Aip 
and (. . ■)mp is quite the correlation function in this model. Corresponding to the usual in 
CFT normalization of fields < ^i,k{x)^i,k{0) >Mf 
1)) we have to take 



\x\ 



Qp 



2{2p-l){p-l] 



p 



p+1 



4^*= (here 4Afe ^ k{kp-2)/{p + 
-, 1/2 



2p 



p + 1 



1 



Thus denoting 



A = 



^p+1 
2^T^Ji'^ 

Qp 



?>p 
P + i 



(4.25) 



(4.26) 
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one can interpret the series (4.19) as the perturbative expansion in the $i3-perturbed M.p, 
i.e., in the RSG model 

Ax = Am,+xJ ^i3ix) (fx (4.27) 

where A_m.^ is the formal action of the CFT }Ap. Under this interpretation the quantity 
(4.19) at 2a/ (3 = 1/p — k becomes the perturbed finite size energy of the primary state 
^i,k+i in the model (4.27). In particular k = and k = 1 correspond to the perturbed 
states = I and $12 respectively. If 2a/P = —1/p + k we can similarly consider first the 
integration over Ui in (4.23) and then treat the remaining integral over Vi as a perturbative 
one in the model (4.27). 

Let us denote F(fc)(p|t) the scaling function corresponding to the finite size energy of 
the RSG state $i,fe+i. If the sin-Gordon coupling /x is related to A as in eq.(4.26) we have 

F^,,){p\t) = FsG{±{l/p-k),p\t) (4.28) 

Comparing with (4.18) we see that the polymer functions (2.7) and (2.8) are (provided 
relation (1.5) between p and n holds) 

F+(n|t)=F(o)(p|t) 
F_(n|0 = F(i)(p|0 

i.e. correspond to the RSG states I and $12 respectively. 



5. Perturbation theory 

Here we present few explicit calculations about the sin-Gordon perturbative series 
(4.19). For the later comparisons it is convenient to use the scaling parameter 

e=(27r)-2(f-i)/(f+i)r7V^p (5.1) 

instead of rj in eq.(4.19) (here Vp is the fx-M ratio (4.8), (4.12)). At p — 2 this notation 
conforms that of eq.(2.17). From (4.19) we have 

F,G{a,p\t) = 5^c2n(«) [(27r)(^-^)/(^+i)rp] (5.2) 

n=0 

At n = 1 the integral (4.23) is carried out explicitly 



C2(a) = 6^ / X (5.3) 

7 



where the notation 7(0;) = r(a;)/r(l — x) is used. 
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In the polymer language eq.(5.2) together with (4.22) and (5.3) predicts the lead- 
ing t — > behavior and the ncxt-to-lcading correction to the polymer scaling function 
F{m, n\t). In particular, expanding around p = 2 and 2a/P = 1/2 we find 

tF(™,n|t) = m(-l + j^« + ...)+n(l-j^« + . ..)+... (5.4) 

in agreement with eqs.(2.19) and (2.30). 

At 2a/(3 — 1/p (this value corresponds to m = n in the polymer language) we can go 
two steps further thanks to the RSG interpretation described in sect. 4. Namely 

C4 (1^ = i) = J^/ {*13(l)'l>13M)A<J«r''/<''+'>d^l. 

^(,,,V(^_1),(-^) (.5) 

2p-2\_^f p 



4(p-l)^(2p-l)^7^— - 7 



p+1 J VP + 1 



and 



{j^D^ 2^ / (^13(l)$13(«l)$13(«2))Mj«l«2r'/(^+')rf'«l(i'«2 

:s(;.:r.(^), -(-^).(^ ) ,.e) 



(p-l)^(2p-l)^7^— - 7 



p+l J \p + l 

At p = 2 + 2n/7r + O(n^) this gives for the RSG scaling function F(o)(p|t) (see eq.(4.28)) 

t^(o,(p|t) =n[-l- + + ■ . .) + 0(n^) (5.7) 

in complete accordance with the expansion (2.31). 

Another special case of the oppositely weighted winding and non-winding polymer 
loops m — —n is related by eq.(4.17) to the sin-Gordon quasimomentum 2a/P = ±(l/p— 1). 
This point again admits the RSG interpretation corresponding now to the finite-size ground 
state in the $12 sector of the $i3-perturbed /Ap. Therefore 

/2q: 1 \ 3 /" (Pu 

C4[ — = --l)= / < $12(0)$13(1)$13(W)$12(00) 

where -C . . . ^m.^ denotes the connected cylinder correlation function in M.^,. Explicitly 
< $12(0)$13(1)$13(«)^12(00) ^M= ai {fi{u)fi{u) - 1) + a2f2{u)f2{u) (5.9) 
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where 



7 



ai 



P 



p + l 



7 



2p- 1 
p + l 



7 



7 



2p 



7 



p — 1 



02 



p+l 



7 



2p 



7 



2-2p 



p + l J \p + I J \p+l 
and /i,2 are in terms of the hypergeometric functions 



7 



p — I 
p + l 



h{u) = (1 - tL^^ J_, 

^ ^ ^ ' \p+i ' p + r p + r J 

f,(u) = _ ^)(2-2p)/(p+l)^ 1^ 

p + l p + l 

The integral (5.8) again can be carried out 

3Y 



u 



2/ 2p-l 



p + l 



2Y 
P 



2p 



7 



P 



p + l J \P + l 
I 



p+l J \p+l 

At p ^ 2 we have therefore 
tF^i){p\t) = n 



p+l 



2(p+l) 



27 



2(p+l) 



18 272(1/3)^ 3274(1/3) 
to be compared with the first terms of eq.(2.32). 



(5.10) 



(5.11) 



+ 0{n^) 



I 
2 

(5.12) 
(5.13) 



6. TBA considerations 

In the previous two sections the sin-Gordon scahng function FgG{a,p\t) has been 
studied in the UV perturbation theory. In principle this approach provides us with a 
systematic expansion of Fsoict, p\t) in powers oft2/(p+i). At the same time the integrability 
of the sin-Gordon model allows the same observable Fsg to be considered in the TBA 
framework. In TBA one does not concern directly the field theory action starting instead 
from the (exactly known in the sin-Gordon case) relativistic factorized scattering theory. 
After some manipulations the cylinder scaling function is related to a system of non-linear 
integral equations (the TBA system), its form depending on the factorized scattering 
amplitudes. 

In the case of the sin-Gordon scattering theory the TBA system (at p > 1) is borrowed 
essentially from the construction by Takahashi and Suzuki [18] for the anisotropic Heisen- 
berg chain. The structure depends drastically on the arithmetic nature of the sin-Gordon 
parameter p. In general it is an infinite system of coupled non-linear integral equations for 
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infinite number of unknown functions. At p rational it can be simplified and reduced to a 
finite system (see [18] for the details). This is why in the TBA studies of the sin-Gordon 
model one tends to choose p rational. It should be mentioned however that this "fractal" 
dependence on p is a problem of the TBA technique itself. The sin-Gordon physics is quite 
continuous in the coupling constant. 

For the time being we are interested in the point p = 2 and its vicinity. At p = 2 the 
TBA system contains only two unknown functions ei (9) and 62 (9) and reads 



t cosh 6* = ei + s * log (1 + 2 cos (47ra//3) e"''^ ^ e"^^^) 
= e2 + s*log(l + e-"i) 



where * denotes the convolution in 9 



s * 



with the kernel 



It happens that at m = 2cos(27rQ;//3) small, ei ~ — logm and 62 ~ m. Denoting 
we find that at m — > 



tcosh^ = e + s * log(l -I- T]^) 
rj = s * 



(6.1) 



/oo 
s{9 - 9')g{9')d9' (6.2) 
-00 



f e'"^' dw _ 1 

i 2 cosh !^ 27r ~ 27r cosh 9 ^^'^^ 
2 

The sin-Gordon scaling function appears as 

FsGia,2\t) = j cosh^log(l + e-^i(^))d^ (6.4) 



(6.5) 



cosh ^e-^(^) ^ +0(m2) (6.6) 
27r 

while e{9) and r](9) solve the following system 



(6.7) 



which is quite the one quoted in sect. 3 eq.(3.21) (with A = l/47r). Note, that eqs.(6.7) 
imply the following functional system for R+{9) — exp{—e{9)) and r]{9) 

R+{9 + m/2)R+{9 - m/2) = 1 + r^\9) 
r]{9 + i7T/2) + r]{9-i7T/2) = R+{9) ^ ' ' 
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Moreover, it can be shown that under appropriate analytic and asymptotic restrictions on 
R+{0) and r]{9) the functional and integral systems (6.8) and (6.7) are equivalent. 

To arrive at eq.(3.23) one needs to shift slightly from the point p = 2. For the reasons 
mentioned above and following ref.[l] we choose the series of rational values 

p = 2 + l/N , iV = 2,3,4,... (6.9) 

Any analytic at p = 2 function of p is unambiguously recovered from the series of its values 
at (6.9). 

The Takahashi-Suzuki TBA system at p = 2 + 1/7V is quoted in [1]. It includes N + 2 
functions ea{0); a = 0,1, . . . , N + 1 and has the form 

t cosh^ = Eq + s * Li 

= Si + s*Lo + sn*L2 + dN * Li 
= £2 + sjv * [La - Li] (6.10) 
3 < a < AT, = ea + SN* [La+i + La-i] 
= sn+i + sn * Ln 

The kernel s{9) is given by eq.(6.3), 

J 2cosh— 27r 
2N 

n{N-l)LO (6.11) 



dN{0) 



e'''^ cosh . , 

2N 

2 cosh — — cosh — — ^" 
2 2N 



while the functions La{9); a — 0,1,2, . . . , N + 1 in eq.(6.10) are read as follows 

Laid) = log (l + e-^'^^^^y 0<a<N 
Ln+i{0) = log (l + ge-^^+i^^)) (l + g-ie-^^+i(^)) 



(6.12) 



Here 

9 = exp(^^(2iV + l)) (6.13) 

introduces the sin-Gordon quasimomentum a. The cylinder scaling function is determined 
by Lo{e) 

FsG{a,2 + l/N\t) = -^ J cosh ^Lo(^)d^-^ tan ^ (6.14) 

For our purposes it is more convenient to consider the functional system for Ya{9) = 
ex.p{-ea{d)), a = 0, 1, 2, . . . , AT + 1, which follows from (6.10) 

io(o+f)io(o-f)=i + i'.W ; 
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2N \ 2N 



2N J J \ \ 2iV 



2N) ^\ 2N J 1 + Yi{9) ' 
2N J " V 2N J 



+ ^) " ^) = + >Ov-i(^)) (1 + glOv+i(^)) (1 + q-'Yj,+,{e)) ; 

(6.15) 

It is easy to verify that this system is "solved" in terms of a single function G{9) (defined 
up to an overall multiplying constant) as 

{q'G2N+i{e) - q-^G.2N-m) {GM - G.m) 



{qG2N+i{e) - q-'Gi{9)) {qG-i{e) - q-^G-2N-i{e)) 

-, , y {Gn+i{0) - G'Ar-l(6')) {G-Nj^i{9) - G-N-l{0)) /g 

{qGM-i{0)-q-'G.M-im{qGN+i{e)-q-^G.N+im ^ 
gGi(^)-g-^G-i(^) 
Gi(^)-G_i(^) 

and for 2 < a < TV 

1 I Y,{e)- {^GN-aie)-q-'G.^+a.2{e)) {qG^.a+2{e)-q-'G_j,+M) ^^^^^ 

{GN-a+2{0) — Gjv-o(^)) {G-N+a{0) — G-N+a-2{0)) 

provided G{d) satisfies 

GQN+2{d) = q-''G{d) (6.18) 
In eqs.(6.16-18) we use the abbreviation 

Gk{e) = G{e + '^^ (6.19) 

Note that due to (6.18) all the functions Ya{9) are m(3 + l/A^)-periodic in 6. 

Turn at the limit N ^ oo and 2a/P = 1/2 — e/2 with e — > 0. Following eqs.(1.5) and 
(4.17) 

m = 7re + O(e^) 



"=277 + ^VAr2 
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1 \ (6.20) 



In this notation q = ) exp(— z7re(Ar+ 1/2)). It is convenient to introduce new function 

m=Ome^{-'-^^^^0) (6.21) 

which enjoys the m, n-independent behavior under the period shift (6.18) 

h{d + iTr(3 + l/N)) ^ -h{e) (6.22) 

The function Yq{9), which is most important for us, now becomes (again hk{9) = h{9 + 
mk/2N)) 

^ (e»V3/,2^+i(^) + e^'5/ii(^))(e-^^/i_i(^) + e-^V3/,_2^_^(^)) ^ 

with d — TTe{N + l/2)/{N + 1/3). At A'" = oo and e = this quantity vanishes as one could 
expect, while h{$) turns to some limiting function g{9) such that 

g{9 + 3i7r) = -g{9) (6.24) 

Moreover, the first order terms in m and n 

Yo{9) = mR+{9) - nR_{9) + higher order terms (6.25) 

are determined by g{9) 

R 2iig{9 + i7r)-g{9-t7r))gi9) 

^ - {g{9 + ZTT) + gi9)) ig{9) + g{9 - z.r)) ^'-"'"^ 
o ... 2i{g{9 + t'K)-g{9-t'K))g\9) 

) (^(^ + Z7r) + <7(0))(<7(0) + <7(^-Z7r)) ^ 

Denoting also 

... . g{9 + in/2)-g{9~in/2) 

^^^^ = ~'WTJ^^J2)TW^J^) ^'-'^^ 

one finds that (6.26a) and (6.27) satisfy the functional system (6.8). Therefore we have to 
identify these functions with the solutions to the TBA system (6.7). 

From (6.27) one can find g{9) in terms of ri{9) 

log g{9) = ~ J tanh(^ - 9') axcia.nr]{9')d9' + const (6.28) 

Relation (3.23) follows. 

The TBA system (6.7) together with (3.23) has been used to compute numerically 

F„(t) = -i- y cosh 9 R+{9)d9 (6.29) 
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and 



for the plots of figs. 5-7. 



Fcit) + 



1 

2^ 



cosh eR-{e)de 



(6.30) 



7. Concluding remarks 

Expansions (2.19) and (2.28) imply that the scaling functions tFn{t) and tFc{t) are 
still real (at least at |t| small enough) after the analytic continuation 



'""IH (7.1) 

(or ^ — > — ^ in eq.(2.1)). As it is argued in [9] such an analytic continuation would 

correspond to the change from the subcritical (dilute) polymer phase to the supercritical 
(dense) scaling polymers. The continued functions Fn(f) = e^^'^ ^ ^ Fn{e^'^^ ^ ^t) and Fc{t) = 
^3i7T / 4 p^^^SiTT / 4^^ are given (inside the convergence region) by the alternated series (2.19) 
and (2.28) (again ^ = t^/S) 



1 4 °° 
tF^(t) = -- + -Y,i-rnUnC 

1 A 



3 3 

(7.2) 



n=l 

with the same Un and Vn as in sect. 2. Outside the convergence region the "dense polymer" 
scaling functions can be found as a particular solution to the continued equation (2.13) 

I ^t^U = -I cosh2U (7.3) 
tdt dt 2 ^ ^ 

fixed by the initial condition at t — > 

t/(t) = -ilog^-ilog« + 0(t^/3) (7.4) 

Now 

t^n{t)=t^ (7.5) 

while Fc{t) is restored from the continuation of eq.(2.23) 

1 1 /•* 

tF^{t) = _ + - / sinh2t/(r) rdr (7.6) 
18 2 Jq 

Contrary to the dilute polymer case where U (t) is always regular at real t > 0, the continued 
function U (t) has an infinite number of singularities on the positive real axis. These appear 
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as double poles of exp(— 2C/(t)) or as simple poles in Fn{t) and Fc{t). In ref.[9] these poles 
were attributed to the level-crossing effect observed in the finite-size dense polymer system. 
Note however that exp(2[/) has double zeroes instead of the poles and, as it is readily 
figured out from eq.(7.3), remains always finite at t > 0. Therefore tF-{t) — tFc{t)—tFn{t) 
(the continued tF-{t) of eq.(2.32)) is positive and non-singular with an infinite sequence 
of critical points. 

Anyhow the analytcally continued Painleve III transcendent U{t) worth more atten- 
tion. Could one invent something as effective as the Predholm or TBA-like representations 
(3.10) and (3.21) for this case? 

A closely related continuation problem arises in connection with the general sin- 
Gordon scaling funciton (4.13). As one concludes from the structure (4.19) the continuation 
r]'^ — > — (or equivalently t e^'^^^+'^^/^t) leads to a real expansion for the function 

tFsG{a,p\t) = tFsG (f^^+^^H) (7.7) 

From the perturbative point of view this alternated series would correspond to a purely 
imaginary coupling constant /i in the sin-Gordon action (4.1). Being defined this way the 
imaginary coupled sin-Gordon model (ISG) is apparently non-unitary and the common field 
theory intuition fails to figure out the structure of its space of states, vacuum etc. However 
one could proceed formally (say perturbatively) arriving at some sensible conclusions. E.g. 
the renormalization group calculations in ISG with 1 formally work and show that at 
least at p large enough ISG is massless and interpolates between two sin-Gordon critical 
points with different values of p. Moreover the local arguments about the sin-Gordon 
integrability are independent on the nature of /U. Therefore the same local higher-spin 
integrals of motion are expected in the ISG as well. In ref.[9] a consistent factorized 
scattering theory for the massless ISG excitations was proposed and supported by the 
calculation of the ground state energy in an external field. There are some difficulties with 
the TEA treatment of this scattering theory and up to now it is not clear if it is complete 
or other (massive) particles are present in the ISG spectrum. 

The interest to ISG is not purely academic. As it was suggested in [9] this field theory 
model is closely related to the dense phase of the 2D polymer problem (in its scaling 
limit). In particular the renormalization group behavior of the ISG effective central charge 
conforms qualitatively that observed in the finite-size dense polymer system [1,9]. 

The TEA approach allows the finite-size dilute polymer scaling function to be evalu- 
ated very carefully. Unfortunately the lack of correct TEA equations impeds the same for 
the dense polymer case. In this connection I'd like to mention a new potential approach to 
the finite-size problem in the sin-Gordon model developped by Destri and deVega (DdV) 
[19]. Contrary to TEA one does not deal with the physical scattering theory, thermal equi- 
librium of physical particles etc., but starts instead with some "constituent particles" and 
their "bare" scattering. After a kind of renormalization DdV arrive at a system of integral 
equations which are free of the "bare" parameters (including instead the "renormalized" 
physical spectra and amplitudes) and resemble strongly the usual TEA equations. It is 
not yet clear if this approach can be generalized to other integrable relativistic models 
or how the Takahashi-Suzuki TEA equations are related to the DdV ones. However the 
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DdV system is verified (numerically) to predict the same finite-size energy as the TBA 
approach does, being at the same time in many respects much more convenient (in par- 
ticular the DdV system is continuous in p, as opposed to the Takahashi- Suzuki TBA). 
Moreover a slight modification of the DdV system looks quite suitable for the analytic 
continuation (7.7) of the sin-Gordon scaling funciton. This modified DdV system will be 
reported elsewhere. 

Finally, the following remark seems in order. The form of potential (3.9) suggests to 
consider a generalization 

2u{e) = . . . + fae-^^ + he-^ + he^ + he^^ + he'^ + . . . + r(^) (7.8) 

It is readily verified that in this case (in the same notations as in eq.(3.8)) 

dWi d'^Wi 



dh dtl ° 

— 



(7.9) 



dh dt\ V dtl 

A natural guess is that in general the function (3.10) satisfies 



-2U^ (7-10) 



dt^ dtl \dti 

i.e. the modified KdV equation. Considering derivatives in t^ etc. we can speculate 
some higher integrable differential equations of the KdV hierarchy, U being related to the 
corresponding r-function. 
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Figure Captions 

Fig.l. Topologically different configurations of two closed polymers. 

Fig. 2. Infinite honeycomb cylinder. 

Fig. 3. Winding and non-winding polymers on a cylinder. 

Fig. 4. Distinct cylinder configurations of two closed polymer loops. 

Fig. 5. Non-contractible one-loop scaling function tFnit). 

Fig.6. Contractible one-loop scaling function tFc{t) + with the bulk term subtracted. 

Fig.7. Scaling functions tF+{t) + t^/A (sohd curve) and + t^/A (dashed one). 
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